home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / slaqsp.z / slaqsp
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSLLLLAAAAQQQQSSSSPPPP((((3333FFFF))))                                                          SSSSLLLLAAAAQQQQSSSSPPPP((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SLAQSP - equilibrate a symmetric matrix A using the scaling factors in
  10.      the vector S
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE SLAQSP( UPLO, N, AP, S, SCOND, AMAX, EQUED )
  14.  
  15.          CHARACTER      EQUED, UPLO
  16.  
  17.          INTEGER        N
  18.  
  19.          REAL           AMAX, SCOND
  20.  
  21.          REAL           AP( * ), S( * )
  22.  
  23. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  24.      SLAQSP equilibrates a symmetric matrix A using the scaling factors in the
  25.      vector S.
  26.  
  27.  
  28. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  29.      UPLO    (input) CHARACTER*1
  30.              Specifies whether the upper or lower triangular part of the
  31.              symmetric matrix A is stored.  = 'U':  Upper triangular
  32.              = 'L':  Lower triangular
  33.  
  34.      N       (input) INTEGER
  35.              The order of the matrix A.  N >= 0.
  36.  
  37.      AP      (input/output) REAL array, dimension (N*(N+1)/2)
  38.              On entry, the upper or lower triangle of the symmetric matrix A,
  39.              packed columnwise in a linear array.  The j-th column of A is
  40.              stored in the array AP as follows:  if UPLO = 'U', AP(i + (j-
  41.              1)*j/2) = A(i,j) for 1<=i<=j; if UPLO = 'L', AP(i + (j-1)*(2n-
  42.              j)/2) = A(i,j) for j<=i<=n.
  43.  
  44.              On exit, the equilibrated matrix:  diag(S) * A * diag(S), in the
  45.              same storage format as A.
  46.  
  47.      S       (input) REAL array, dimension (N)
  48.              The scale factors for A.
  49.  
  50.      SCOND   (input) REAL
  51.              Ratio of the smallest S(i) to the largest S(i).
  52.  
  53.      AMAX    (input) REAL
  54.              Absolute value of largest matrix entry.
  55.  
  56.      EQUED   (output) CHARACTER*1
  57.              Specifies whether or not equilibration was done.  = 'N':  No
  58.              equilibration.
  59.              = 'Y':  Equilibration was done, i.e., A has been replaced by
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSLLLLAAAAQQQQSSSSPPPP((((3333FFFF))))                                                          SSSSLLLLAAAAQQQQSSSSPPPP((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              diag(S) * A * diag(S).
  75.  
  76. PPPPAAAARRRRAAAAMMMMEEEETTTTEEEERRRRSSSS
  77.      THRESH is a threshold value used to decide if scaling should be done
  78.      based on the ratio of the scaling factors.  If SCOND < THRESH, scaling is
  79.      done.
  80.  
  81.      LARGE and SMALL are threshold values used to decide if scaling should be
  82.      done based on the absolute size of the largest matrix element.  If AMAX >
  83.      LARGE or AMAX < SMALL, scaling is done.
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.